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ABSTRACT 

The nuclei of a wide class of active galaxies emit broad emission lines with widths 
at half maximum (FWHM) in the range 10 3 — 10 4 km s _1 . This spread of widths 
is not solely a consequence of the range of the luminosities of these sources since a 
plot of width versus luminosity shows a large scatter. We propose that the broad line 
emission region (BLR) is axially symmetric and that this scatter in line width arises 
from an additional dependence on the angle of the line of sight to the axis of the 
emission region. Such a relation is natural in unified models of active nuclei which link 
a variety of observed properties to viewing angle. Adopting a simple form for the line 
width as a function of luminosity and angle, and convolving this with the observed 
luminosity function, allows us to predict a line width distribution consistent with 
the available data. Furthermore, we use the relation between the equivalent width of a 
line and the luminosity in the continuum (the 'Baldwin Effect') to predict an observed 
correlation between line width and equivalent width. The scatter on this correlation 
is again provided by angular dependence. The results have applications as diagnostics 
of models of the broad line emission region and in cosmology. 



Key words: 

theory 



galaxies: active: - galaxies: Seyfert - quasars: emission lines - cosmology: 



1 INTRODUCTION 

In unified models of active galactic nuclei with spherically 
symmetric BLR the width distribution of the broad emis- 
sion lines cannot be accounted for by luminosity dependence 
alone. Plots of line width versus continuum luminosit y have 
a large sca t ter and show no significa nt correlation (Wills 
et al. 1993; Puchnarewicz et al. 1997). There is, however, 



(iii) Several axisymmetric disc- wind models, su ch as those 
of Cassidy & Raine ( |l996[ ), Chaing fc M urray ( |l996[ ) and 
Emmering, Blandford & Shlosman ( 1992 ) have been pro- 
posed and mode ls of this type are gaining support from evi- 
dence for winds (Hamann 1997). These will naturally predict 
some viewing angle dependence of line width. 



growing evidence that the broad line region (BLR) is not 
spherical, but axisymmetric. 



(i) Observed samples of AGN (Wills and Browne 1986 



Wills et al. 199E ; Brotherton et al. 1994| ) suggest relations 
between line widths an d R, the ratio of core to lo be dom- 
inance. Other samples ( Puchnarewicz et al. 199?] ) find re- 
lations between line widths and a ox , the continuum slope 
parameter from the optical to X-ray bands. Both of these 
parameters have some viewing angle dependence. 

(ii) The continu um and line light curve s of some active 
nuclei, eg 3C390.3 (Wamsteker et al. 1997), are most natu- 



rally interpreted in terms of a disc-like line emission region. 
It has also been suggested that some double peaked line 
profiles arise from discs (eg Arp 102b) although the inter- 
pretation in these cases is not so clear when time variability 
is taken into account. 



It should be noted that Osterbrock (1977) showed that a 
deficit of systems with narrower lines ruled out pure disc 
models, but such objections do not necessarily apply to ax- 
isymmetric models in general. 

In this paper we shall adopt a simple dependence of line 
width on both viewing angle and luminosity. Then: 

(i) We obtain a reasonable fit to the distribution of line 
widths. 

(ii) Given the Balwin relation between line and contin- 
uum luminosity we predict a relation between line width 
and equivalent width compatible with the observed trend. 
The scatter on this relation is attributed to angular depen- 
dence. 

(iii) We discuss how the width distribution can be used 
to test models of the BLR. 

(iv) If the BLR is indeed axisymmetric, we show how the 
line width distribution can be used, in principle, to deter- 
mine cosmological parameters. 
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2 MODELLING THE LINE WIDTH 
DISTRIBUTION 

For an axisymmetric BLR we can expand the line width as 
a function of angle i between the axis and the line of sight 
in a Fourier series 



v(i) = ao + a n f 



with the cosine terms set to zero for symmetry about « = ■§ • 
However the line widths do not appear to be a function of 
angle only as it is difficult to obtain the range of widths re- 
quired in an uncontrived w ay. Also, observat ions of variable 
sources such as NGC 4151 (Ulrich et al. 1991) show that line 
width is luminosity dependent. However, there is no reason 
to assume that the variation in line width in one source 
should be the same as that between sources at different lu- 
minosities. In fact in the case of line equ ivalent widths this 
is definitely not the case (Peterson 1997). Therefore we in- 
troduce the luminosity dependence as another independent 
variable. Thus, we consider v(i, L) to depend on angle and 
luminosity with a n — a n (L), which we expand as 



a n (L) = L' 



m=0 



% L m for n > 0. 



Note that in principle a could depend on n, but we shall take 
it to be a constant for a given line. To lowest non-trivial order 
we terminate the series at aio and write aoo = a, aio = b. 
Thus, we arrive at our proposed model: 



v(i,L) 



(a + b sin i)L° 



More precisely, we suggest that line width at half maximum 
(FWHM), v km s _1 , is related to the B band luminosity, in 
units of 10 44 erg s _1 , and viewing angle as follows 



v(i, L44) = (a + fesini)Z/44- 



(1) 



To produce the line width distribution we combine the 
above relation with statistical information about the number 
of systems at each viewing angle and each luminosity. The 
number of systems, T, at each viewing angle satisfies 

dT(i) = sini di, < i < u < 90°. 

The angle i* takes account of the obscuring material in the 
unified picture that prevents direct sight of the BLR at large 
inclinations. This value will vary between systems and may 
be correlated with luminosity. Realistic unified models prob- 
ably have i* — 2*^44) but we shall take i» = constant, since 
this is all that can be justified by current data. 

The number of systems at each luminosity is found from 
the luminosity function $(L). We use a broken power law 
fit to the lumin osity function given in Boyle, Shanks and 
Peterson (1988) to give 



$(L 4 



OC ^44^ 



p = 3.85 for L44 > 2.0 
p = 1.27 for L44 < 2.0 



The number of systems at each width, N(v), is given by 

poo pi* 

N(v)= / 5(«-w(i,L 44 ))dT(i)$(L44)dL 4 4 
Jo Jo 

wher the 5 - function restricts the integral over luminosity 
and angle to systems with width v. Thus evaluating the 
angular integral, we obtain 



N ( v ) = I T^7T*( i 44)dL 44 . 
'0 



I dv I 

di I 



(2) 



To compute this integral numerically we restrict L44 to the 
range 0.01 < L 44 < 13.0. 



2.1 Equivalent widths 

To complete the picture we introduce the parameter j3 for 
the Baldwin relation between the equivalent width (EW) of 
a line and the blue band continuum luminosity of the active 
nucleus: 

EW oc Lf 4 . 

From the observational data sets we can obtain a relation 
between the EW and FWHM of a line of the form 

EW oc v 1 . 



Using our previous relation, v cc L44, we see 
EW ocv l3/a . 



(3) 



with a scatter deriving from the angular dependence. Thus, 
any two of a, (3 or 7 should determine the third. In particu- 
lar, given the Baldwin relation and a fit to N(v), determining 
a, we predict the slope of the correlations between FWHM 
and EW (with a scatter from angular dependence). 



3 RESULTS 

3.1 Line width distribution 

We have investigated the Civ A1549, Mg n A2 798 and 
H/3 A4861 li nes us ing data from W ills et al. ((t99§), Broth- 
fL994p, Wilkes (|l987[) and Puchnarewicz et al. 



crton et 



; al. (|1994 
Fig.pT shov 



(1997). Fig.|l| shows the observed width distribution in Civ 
for three samples with a theoretical curve drawn for com- 
parison using i* = 60° and the other parameters as shown 
in table 1. This value for i* is chosen to give a reasonable 
fit for the three emission lines simultaneously. This value is 
larger than expected for the standard unified model of AGN 
but co nsistent with the unified model of Cassidy and Raine 
( [l996| ). The position of the peak of the distribution is depen- 
dent mostly on a, b and i, with a determining the spread of 
the distribution. Setting a = gives a pure disc model and 
increases the predicted number of narrower lined systems. 
Setting 6 = gives a spherically symmetric BLR model and 
the distribution is a narrow peak centred around a (assum- 
ing a is small). The Wilkes data, taken from the Parkes flat 
spectrum survey, is biased towards face on objects, which 
this model predicts would give a peak in the distribution at 
lower velocity. The data in Wills et al. and Brotherton et al. 
is biased towards higher luminosity objects which may ex- 
plain the lack of narrower lined systems in the lower figure. 
Fig.p shows similar fits to Mgn and H/3. 

3.2 FWHM - EW 

Fig.[] shows the corresponding EW - FWHM relations 
obtained from the model. The slopes, from table [l], are 
/3/a = -0.85,1.1,0.89 for Civ, Mgn and U/3 respectively 
compared with — 1±0.25 for Civ (from inspection of Wills 
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Table 1. Model parameters used to fit the observed samples for 
3 emission lines. Value of f3 taken from Peterson (1997) for CIV 
and from Puchnarewicz, E.M. (private communication) for Mgii 
and H/3 



Line 


a 


b 


a 


P 


Civ 


2500 


9500 


0.15 


-0.17 


H/3 


2500 


8000 


0.30 


0.4 


Mgn 


500 


2500 


-0.35 


-0.37 



viewing angle increases. Whilst this is in the opposite sense 
to the other models we can still obtain the line width distri- 
bution using equation (jj]) with b < and so cannot rule out 
the magnetic winds model. However where correlations exist 
between FWHM and R or a ox , they suggest that v increases 
with i, hence b > 0. It may be possible to discriminate be- 
tween the models by consideration of the predicted value of 
a. 



et al. fll993|) ), 1.6±0.4 for Mgn and 0.6± 0.4 for H/3 
(Puchnarewicz, E.M., private communication). The sharp 
edge to the data in these figures will be disrupted by time 
variability in the source luminosity since the Baldwin rela- 
tion in a single variable source does not h ave the same slo pe 
as the relation between different sources (Peterson 1997). 



4 DISCUSSION 

4.1 Mgii an exception? 

Whilst it is perhaps unexpected for the Mgn emission line 
distribution to be fitted with a < 0, this is necessary for con- 
sistency with the observed Baldwin effect (/3 < 0) amd the 
positive correlation between EW and FWHM {ft /a > 0). Al- 
though a positive a gives a similar distribution curve, a < 
is physically more reasonable as it corresponds to a decrease 
in Mg + ions with increasing ionization parameter. Systems 
which have comparably wide Mg n and C iv lines are rare 
and seem to occur at low luminosity (e.g. NGC 4151). More 
commonly the widths vu ne satisfy fJMgii <C ^civ- Presum- 
ably at high luminosity the low ionization emission lines are 
absent from the inner, faster moving, BLR material. Cer- 
tainly at higher luminosity (e.g. NGC 5548) reverberation 
mapp ing shows Mgn t o come predominantly from larger 
radii flClavel et al. 199l|). 



4.4 Application to cosmology 

Since this analysis enables us to remove statistically the 
scatter in the luminosity-FWHM relation for AGN it has 
obvious applications to cosmology once large data sets of 
FWHM become available. Note that the observed distribu- 
tion of line widths is itself independent of any cosmological 
model. If other methods can be used to determine the line of 
sight to individual objects these too could then be used for 
cosmology. In further work we expect to include a redshif t 
dependence in the luminosity function $(L, z) (Pei 199E) 
with L — L(Hq, go), to provide a new method of determin- 
ing Ho and q . 



5 CONCLUSIONS 

We conclude that the simple picture we have presented here 
accounts for the scatter in FWHM versus luminosity, ac- 
counts for the distribution of FWHM, and relates the trend 
in the EW with FWHM to the Baldwin relation. This may 
be useful as a diagnostic tool in discriminating between disc- 
wind models. The analysis has applications as a cosmological 
tool particularly as measurement of line widths is indepen- 
dent of any cosmological model. 



4.2 Narrow lines 

Narrow lined systems (FWHM< 500 km s _1 ) appear to be 
under-represented in the samples. This may be a resolution 
or a selection effect. However a pure disc model (a — 0) 
would pre dict l arger numbers of such systems. This led Os- 
terbrock ( 1977 ) to rule out disc models for the BLR. How- 
ever, it is as yet unclear where the Narrow Line Seyfert 1 
class of o bjects fits into the scheme. Rec ent analysis of these 
systems (Rodriguez-Pascual et al. 1997) shows the existence 
of a broad component to the narrow lines suggesting that the 
BLR is not obscured in these sources, which may therefore 
provide a class of "broad lined" objects with small FWHM. 



4.3 Comparison with theoretical axisymmetric 
models 



Current axisy mmetric models of the B LR include the dual 
winds picture ( Cassid y and Raine 1996 ) , the disc wind model 
( Chaing and Murrayl996| ) and the magnetic wind model 
(Emmering et al. 1992). The dual winds model and the disc 



wind model both have an angular dependence of FWHM on 
viewing angle which is consistent with our model (Figjij). 
Examination of figure 6 in Emmering et al. (1992) suggests 
that in the magnetic wind model the FWHM reduces as 
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Figure 1. Model predictions for CIV overlaid with observed data by Wilkes (1987) (top) and Wills et al. (1993) (dashed) and Brotherton 
ct al. (1994) (dot-dashed). 
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Figure 2. Model predictions for H/3 (top) and Mgll (bottom) overlaid with observed data from Puchnarewicz et al. (1997). 
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Figure 3. FWHM - EW relation as predicted by this model for C IV. Notice the large scatter produced by the angular dependence 
v. The sharp edge is expected to be disrupted by source variability. 
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Figure 4. FWHM - viewing angle relation predicted by the dual winds model (top) and the disc wind model (bottom). Both models 
are fitted well by a sine curve. 
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